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1 Introduction 

Observing astropliysical neutrinos opens a new window on the Universe. The claim is already 
justified by measurement of solar neutrinos and by Supernova 1987A. Several new detectors 
are in operation or under construction to study in detail the energy range from a few MeV 
to tens of GeV. The primary focus of a kilometer-scale detector will be to observe neutrinos 
at higher energy. 

After reviewing the energetics of extragalactic cosmic rays, which are closely related to 
energetic neutrinos in most scenarios, I list some possible sources that have been suggested 
and describe their neutrino signatures. Next I discuss the nature of the neutrino signals and 
how they can be detected. The summary section includes some remarks on the prospects 
for actually finding signals of high energy astrophysical neutrinos. I also comment on the 
backgrounds of atmospheric neutrinos, which may also serve as a calibration source. 

Apart from a brief mention of diffuse galactic neutrinos, I focus in this talk on neutrinos 
associated with the highest energy cosmic rays, which are assumed to be of extragalactic 
origin. I also do not address here the important subject of neutrinos from decay of neutrali- 
nos, which was covered in David Schramm's talk. (See also Ref. for a recent treatment 
of the subject.) 



2 Extragalactic cosmic radiation 

High energy neutrinos require the presence of energetic hadrons. It is therefore of interest to 
use the observed energy density in cosmic rays to estimate the level of neutrinos that might be 

*Talk given at the OECD Megascience Forum Workshop, Taormina, Italy, 22/23 May 1997. Work 
supported in part by the U.S. Department of Energy. 
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expected. I discuss the extragalactic component of the cosmic radiation and corresponding 
extragalactic sources of neutrinos. An analogous argument can be made for sources inside 
the galaxy 0. 

Fig. shows the cosmic-ray spectrum with an assumed extragalactic component fitted 
above the "ankle" and a GZK cutoff around 5x 10^^ eV. The energy content of this component 
is ~ 3 X 10~^^ erg cm~^, assuming an extrapolation to low energy with an E'"^ differential 
spectrum. (The total energy content would be higher if the spectrum of the extragalactic 
component below the ankle is steeper.) The power required to generate this energy density 
in 10"^° yrs is ~ 3 x 10^'' erg/s per (Mpc)^. This works out to 

• ~ 3 X 10^^ erg/s per galaxy, 

• ~ 3 X 10^^ erg/s per cluster of galaxies, 

• ~ 2 X 10'^^ erg/s per active galaxy, or 

• ~ 2 X 10^^ erg per cosmological Gamma Ray Burst (GRB). 

The coincidence between these numbers and the observed energy output of these objects in 
other wavelengths, especially in the last two cases, is the basis for their use as models of 
origin of the highest energy cosmic rays. 

Suppose now that there is the same energy density in neutrinos, p^; ~ 3 x 10^^^ erg/cm"^ 
with a spectrum 

= AE-^ cm-^s-^sr-^ (1) 

dEy 

that continues to a maximum neutrino energy, E^ax- Then 

J AE-^dE, = 1^^^ (2) 

normalizes the (diffuse) neutrino flux. The resulting signal depends on the neutrino spec- 
trum. For 7 = 1 and E^ax ~ 10* GeV the rate of neutrino-induced muons is about 50 
per km^ per year and depends only weakly on Emax- For 7 = up to the same Emax the 
corresponding number is 20 events per year, with the rate nearly inversely proportional to 

Emax • 

3 Possible sources and their neutrino signatures 

Is there a physical realization of a source that would produce approximately equal energy 
densities in cosmic-rays and neutrinos? Consider a generic source in which protons and 
electrons are accelerated to high energy at shocks. The acceleration process is fast but 
the charged particles remain trapped in the diffuse, turbulent magnetized plasma until they 
radiate photons, the electrons by bremsstrahlung and the protons by photoproduction. Neu- 
trons from p + 7 — >■ 77, -|- 7r° escape to become cosmic-rays. Neutrinos from charged pions will 

iProm Ref. D with the addition of data from Tibet [||. 
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Table 1: Diffuse neutrino rates 



Source 


u jj, 

per km^ yr. 


(TeV) 
(median) 


Z/e + Ue 

(total rate per km^) 


z/gC only 


Atmospheric 


5 X 10^ 


0.1 






Neutralinos 


< 10^ 


< 0.1 






Prompt /i 


400 — 10^ 


0.4 






Galactic center 


300-1000 


0.4 






Galaxy clusters 


< 5000 


1.8 


< 3000 


< 50 


GRB sources 


200 


7 


40 


10 


AGN jets 


500 


15 


70 


200 


AGN cores 


< 3 X 10^ 


2 


< 5000 


< 25 


GZK 


< 1 


500 






Top. Defects 


< 1 


5000 







have an energy comparable to the energy in escaping cosmic-rays. The neutrino/7-ray ratio, 
among other things, reflects the ratio of accelerated electrons to protons in the source. This 
is a crude generic description of some models of AGN jets [^, 

Many of the suggested sources of extragalactic cosmic rays including AGN 0, clusters 
of galaxies j^, ^, |lO], 11], cosmological gamma-ray bursts |12| and radiation from topological 
defects 
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also imply a correlated flux of high energy neutrinos. In addition, there will 
be neutrinos at some level from cosmic-ray interactions on the 2.7 °K background radiation 

mi- 

In Fig. 2 and Table 1 I summarize the event rates in a kilometer-cubed detector from 
various sources of high-energy neutrinos. The figure displays the rates for neutrino-induced 
muons and the table includes events generated by electron-neutrinos as well, including sep- 
arately the rate of events at the Glashow resonance. The diffuse fluxes referred to here are 
obtained by integrating over all sources in the Universe for each model (e.g. all blazars or 
all GRB, etc.). A discussion of the relation between signals from individual point sources 
and the corresponding diffuse fluxes will be made elsewhere |]15[. Generally, the signal rate 
from a nearby point source will be significantly lower than the corresponding diffuse signal, 
as illustrated by the example of 3C273 in Ref. H]. 



4 Neutrino signals 

It is clear from Table 1 and Fig. 2 that some discrimination on neutrino energy will be needed 
to see the signals of high-energy neutrino sources above the local atmospheric backgrounds. 
The rate of neutrino-induced muons is given by 

dRate . dN^ „ ,^ ^ , , , 

Zh^ = " d^^ " 



3 



where 



rE[u] da (E ) 

= NA^^^dE, (4) 

~ NAcy,{E,)Range[{l- {y))E,], (5) 

where Na is Avogadro's number and y is the fraction of energy to hadrons in a charged cur- 
rent interaction of a (or v^) . The last hne of this equation allows one to make reasonably 
accurate estimates easily which display the dependence of the signal on neutrino energy, as 
in Fig. 2. To display the dependence on muon energy it is straightforward to include the 
convolution of the differential cross section with the neutrino range as obtained, for example, 
from Ref. 



For results here I have used the cross sections tabulated by Gandhi et al. IjT^- The 
corresponding interaction length of neutrinos is shown in Fig. 3. Note that the Earth, 
which is 1.1 X 10^ km.w.e. along a diameter, partially obscures upward neutrinos with 
Ej^ > 100 TeV. Fig. 4 is an example of the convolution of a predicted diffuse neutrino flux 



12 | with the projected area (assumed to be 1 km ) times P,^. The bottom panel shows the 
differential rate integrated over the entire upward hemisphere with the effect of absorption by 
the Earth represented by the shading. Note that this effect leads to an angular dependence 
for neutrinos with E > 100 TeV because the vertically upward neutrinos are absorbed. 

For illustrative purposes, I use a simplified treatment of the range that displays the 
essential physical effects of muon propagation. The energy-loss rate is 

with ecriticai = «C ~ 500 GeV. The first term in Eq. |^, which varies slowly with energy, rep- 
resents ionization loss, while the last term represents radiative losses due to pair production, 
bremsstrahlung and hadronic interactions. At the critical energy e the two energy loss terms 
are equal. The solution for the avearage energy is 

E{X) ~ (Eo + e)exv[-X/i] - e, (7) 

and the range R (neglecting fiuctuatons) is obtained by setting E{R) = 0: 

R ^ e ln(l + Eo/e), (8) 



with ^ ^ 2 km water equivalent. (See Ref. [16] for a nice discussion of the effect of 
fluctuations, the effect of which depends on the neutrino spectrum. Formulas and tables 
for muon energy loss in various media are given in Ref. [0.) 

The muon range is shown in Fig. 5 as a function of its initial energy. It is interesting 
that for Efj^fi > 0.5 TeV R> 1 km. Thus for TeV (and higher energy) neutrinos the effective 
volume for neutrino-induced muons is larger that the physical volume of the detector even 
for a kilometer-cubed detector. This is the reason that this mode of detection is of such 
great interest. 
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A charged-current interaction of a (or z/^) can occur anywhere between the full range 
away from the detector to just outside (or perhaps inside) the detector volume. Thus the 
distance a detected muon has traveled before reaching the detector has a median 

Ri ~ ^ein(l + Eo/e), (9) 

with a flat distribution. The corresponding distribution of energies at the detector is therefore 
very broad with a median 

Ef,{median) ^ ey'(l + Eo/e) - e. (10) 

The right-hand curve in Fig. 5 represents this relation. For example, the median energy at 
the detector of a muon that starts with 4 TeV is 1 TeV, while for 200 TeV it is 10 TeV. For 
i?^ ^ e E ^{median) ^ y/eEo. 

As a specific example, I consider the blazar model of diffuse fluxes of Protheroe [0. 
Fig. 6 shows the differential rate on linear (upper panel) and logarithmic (lower panel) 
scales. Absorption in the Earth has been included. The thin solid line is the background 
of atmospheric neutrinos exclusive of prompt neutrinos. The dashed lines give minimal and 
maximal estimates of the prompt neutrino background, which I will discuss further below. 
A similar model by Mannheim P] produces a similar, but somewhat lower, signal. 

A typical neutrino-induced muon from this spectrum will have initial energy of ~ 1000 TeV, 
for which the median energy at the detector is ~ 20 TeV. What will such a muon look like 
as it propagates through the detector? The answer is shown in Figs. 7 and 8, made from 
calculations by Todor Stanev. They show the propagation history of two individual muons 
through 1 km of water. In each case the upper panel is on a linear scale and the lower is 
the same event on a logarithmic scale. What is plotted is the total number of electrons and 
positrons along the track. Each burst generates an air shower with typical length in water 
of 5 m. Since for each subshower 

E shower ~ 1 GeV X Size at maximum (11) 

for showers in this energy range, the vertical axis can also be read as number of GeV deposited 
per burst. 

The examples shown are for E^ = 10 TeV and Efj_ = 100 TeV at the detector. The 
radiative processes approximately scale with energy in this region, so we can summarize the 
situation by noting that a muon with E^ ^TeV will typically generate a couple of bursts 
per kilometer in which it radiates about 10 per cent of its energy. This is the characteristic 
signature that a km^ detector should be sensitive to. 



5 Discussion 



5.1 Some other sources 



The original estimate of diffuse neutrinos from cores of radio quiet AGN [|r^, ^] is still 
of interest. Though the highest predictions of Szabo and Protheroe are above the Frejus 
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limit 1^ and the possibility of neutrinos from radio quiet AGN has been overshadowed by 
the spectacular results from Whipple, HEGRA etc. Q on high energy gamma-rays from 
blazars, the model itself could still be valid. The upper limit allowed by Frejus is indicated 
by the upper dotted line in Fig. 2. 

At a lower level, in addition to the diffuse fluxes from AGN jets, there could be production 



of neutrinos by cosmic rays interacting with gas in clusters of galaxies |TD[. The upper 
dashed curve in Fig. 2 is the signal from maximum allowed flux of Ref. if the same 
source accounts for the diffuse extragalactic 7-ray fluxes above 100 MeV. Estimates based 
on cosmic-ray power of current normal galaxies and AGN are considerably lower than this 
maximum. Achieving this upper limit would require a significant "bright phase" of cosmic- 
ray production at early epochs. 

GRB sources could also produce an interesting flux of diffuse neutrinos if they are the 
sources of extragalactic cosmic rays [1^, as shown by the lowest dotted curve in Fig. 2 and 
also in Fig. 4. 

Another possibility is our own galactic center. This is of special interest because of the 
recent Egret observation that the spectrum from the central region of the galaxy is harder 



than would be expected if the 7-rays were produced from cosmic-rays with the differential 
spectrum of -2.7 observed locally. The lowest dashed curve in Fig. 2 is normalized to the 
Egret data at ^ 2 GeV extrapolated to higher energy with a differential index of 2.5 and 
assuming a similar spectrum of neutrinos (from decay of charged pions). A more quantitative 
fit ||2^ to the whole spectrum, including bremmstrahlung photons as well as photons from 



decay of 7r°s, gives an even harder spectrum at high energy and a correspondingly higher 
neutrino signal. 

5.2 Atmospheric neutrinos: a calibration source 

The idea here is to use the atmospheric neutrinos and their known angular dependence 



27|, E3 to calibrate the neutrino detector. See Fig. 9 for examples of the angular dependence. 



Since high energy diffuse sources may also have a similar angular dependence (depending 
on their energy spectra) due to absorption in the earth, some selection on energy will be 
important. There is a useful discussion of the interplay between energy threshold and angular 
dependence of high energy signals and lower energy background in Ref. which also 



discusses the contribution of neutrinos from decay of charmed particles to the atmospheric 
background. 

5.3 Charm 

The limiting factor, both for using the atmospheric neutrinos for calibration and in the search 
for diffuse sources is the uncertainty in the component of prompt neutrinos. The paper of 



Thunman et al. [30 gives a thorough review of prompt neutrinos, which come primarily 



from decay of charmed hadrons.| The work of Volkova et al. |^ gives another view (larger 



^Because of the steep rise of the structure funcions at smaU a;, heavy flavors are always produced pre- 
dominantly at threshold in hadronic collisions, so production of beauty will always be suppressed relative to 
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cross section) of prompt neutrinos. The lower dashed hne in Fig. 6 gives the contribution of 
prompt neutrinos from charm as calculated in Ref. assuming only gluon-gluon fusion to 



be effective in producing charm. The upper dashed curve indicates the level of charm that 
has been suggested by others if there is significant charm content in the proton |]31]]. (The 
same two curves are shown in Fig. 2 as a pair of light, solid curves.) 

5.4 The highest energy neutrinos 

Neutrinos associated with the GZK cutoff and especially neutrinos from decay of topological 
defects have extremely high energy, though the fluxes are low. Detectors with effective 
volumes significantly larger than a kilometer cubed would be desirable, e.g. radio detection 



32| , |33| , p4| or air shower detectors with gigantic acceptance ^ . 
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FIGURES 



Fig. 1. The high energy cosmic-ray spectrum. The assumed extragalactic component is 
indicated by the sohd hne above 10^^ eV. 

Fig. 2. Differential rate of //-induced, upward muons per logarithmic interval of neutrino 
energy (per km^-yr). The fluxes are integrated over the full 2it steradians from below the 
horizon, including the effect of neutrino absorption in the Earth. 
Heavy solid line: atmospheric neutrinos from vr- and K-decay; 

Thin solid lines (< 10^ GeV): "miminum" |30] and "maximum" [31| estimates of prompt 
neutrinos; 



Upper dotted line: u from AGN cores-maximum allowed by Frejus limit [^] in calculation 
of 



Lower dotted line: v from AGN cores without confinement [|T9| ; 
Thin solid hne (peak at 10^ GeV) AGN jets from Ref. ||]; 

Thin dashed hne: maximum allowed from clusters of galaxies if this is source of diffuse 
(extragalactic) 7-ray background |1TT|| ; 



Heavy dashed line: Galactic plane in the direction of the galactic center (0.73 sr) ||25 
Heavy dotted line (peak at 10^ GeV): GRB model of Ref. P 



Upper thin sohd line (> 10^ GeV): neutrinos from topological defects |T3| ; 

Lower thin sohd hne (> 10^ GeV): example of neutrinos from interactions of cosmic rays 

with the microwave background radiation (maximum red shift = 2; bright-phase index = 2) 



14 . 



Fig. 3. Neutrino interaction length vs. energy (from Ref. jl 

Fig. 4. Convolution of the predicted GRB flux of z/^ -|- z/^ ||T^ with P^x detector area 
(upper panel, arbitrary units) and the corresponding response curve (lower panel), which 
shows the relative contribution of neutrinos of different energies to the predicted signal of 
upward muons. The dashed curve in the lower panel shows the signal when absorption of 
high-energy neutrinos in the Earth is neglected, while the solid curve includes the effect of 
absorption. 

Fig. 5. Muon range vs. energy (left solid curve). The right curve illustrates the relation 
between muon energy at production and at the detector (see text). 

Fig. 6. Differential rate of neutrino-induced muons per logarithmic interval of parent neu- 
trino energy for the AGN jet calculation of Ref. 0. The thin solid line is the atmospheric 
background from decay of pions and kaons and the dashed lines are minimal |^ and max- 



imal |31[ estimates of prompt neutrinos. The same information is shown in both panels 
(upper with linear scale and lower with logarithmic). 

Fig. 7. History of a single 10 TeV muon propagating through one kilometer of water. Scale 
can be read as GeV per burst (see text). Upper and lower panels show the history of the 
same muon on linear and logarithmic scales respectively. 
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Fig. 8. Same as Fig. 6 for a 100 TeV muon. 

Fig. 9. Angular dependence of upward, z/^-induced muons. The same information is shown 
on linear and logarithmic scales in the upper and lower panels. In each case the upper curve 
is for > 1 TeV and the lower for > 3 TeV. The scale is events per year per bin of 
Acos^ = 0.1 in a projected area of 2 x 10^ m^. 
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3. Neutrino interaction length (knn.w.e.) 

lO^P 
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10 10 10 10 10 10 10 10 10 10 

Energy (GeV) 



2. Rates of i/ — induced fj, per sq. km per year 




E, (GeV) 



4. GRB flux times acceptance 




4 5 6 7 8' 

10 10 10 10 10 10 




4 5 6 7 8' 

10 10 10 10 10 10 



E, (GeV) 



7. SIZE for 10 TeV muon 
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8. SIZE for 1 00 TeV muon 
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